Determination of a flow generating a neutral magnetic mode 
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The problem of reconstruction of a flow of conducting incompressible fluid generating 
a given magnetic mode is considered. We use the magnetic induction equation to derive 
ordinary differential equations along the magnetic fleld lines, which give an opportunity 
to determine the generating flow, if additional data is provided on a two-dimensional 
manifold transversal to magnetic field lines, and show that an arbitrary solenoidal vector 
field can not be a neutral magnetic mode sustained by any flow of conducting fluid. 

1. Introduction 

According to the modern scientiflc paradigm, magnetic flelds of astrophysical objects, 
ranging from planets to galaxies, are often sustained by conducting fluid flows, driven by 
convection in the melted medium in their interiors [12, 11, 20, 16, 19]. These processes 
are governed by the Navier-Stokes and magnetic induction equations (supplemented by 
other equations, such as heat equation and rheology relations, as appropriate). However, 
it is diflicult to study them numerically because of the extreme parameter values involved, 
which require prohibitively high resolution of simulations. Thus, application of analytical 
or semi-analytical methods to the study of astrophysical dynamos appears unavoidable. 
In the present paper we suggest an approach, in principle enabling one "to separate" the 
two fundamental equations; hopefully, this can be useful for investigation of asymptotics 
of astrophysical dynamos. 

Usually the magnetic induction equation is employed for investigation of the evolution 
of magnetic fleld for a given flow of incompressible conducting fluid (which is predeflned 
in kinematic dynamo problems, or supposed to evolve simultaneously when nonlinear 
dynamos are studied). We consider here an inverse problem, investigating which conse- 
quences existence of a neutral magnetic mode bears upon the generating flow. We show 
how the flow can be reconstructed uniquely up to the data which must be provided on 
a two-dimensional manifold transversal to magnetic fleld lines. We demonstrate that an 
arbitrary solenoidal vector fleld can not be a magnetic mode sustained be any flow of in- 
compressible fluid, unless the fleld satisfles a consistency equation in the fluid volume. We 
hope that such an analysis may be useful, in particular, for examination of asymptotical 
properties of various steady magnetohydro dynamic systems and their stability. 

In recent numerical studies of nonlinear magnetic dynamos acting in plasma [8] and 
fluid [13, 5, 6] flows with a prescribed forcing, as well as in thermal convection in a 
horizontal layer of conducting fluid rotating about a vertical axis [25] or in the absence of 
rotation [14, 15], it was discovered that temporal evolution can result in emergence of a 
steady state with a non- vanishing magnetic fleld. Magnetostatic equilibria in ideal plasma 
were discussed in [1]. A magnetic fleld of a steady conflguration is a neutral magnetic 
mode, i.e., a vector fleld belonging to the kernel of the magnetic induction operator. 
Neutral magnetic modes play an important role in large-scale dynamos [21]-[24]. 



We therefore focus on neutral magnetic modes in our analysis. However, a straight- 
forward modification of our approach can be applied for reconstruction of flows for eigen- 
functions of the magnetic induction operator associated with any given eigenvalue, or for 
arbitrary evolving magnetic fields. 

2. Reconstruction of flows 

Consider the magnetic induction equation 

dth ^ riV% + V X (uxh). (1) 

Magnetic field is solenoidal: 

V • b = 0. (2) 

In a steady state magnetic field is a neutral mode of the magnetic induction operator. 
For a given fiow b and molecular difi^usivity rj the operator is elliptic. If magnetic field 
generation in a bounded volume of fiuid is considered and regular boundary conditions 
for magnetic field are imposed, it has a discrete spectrum, with the eigenvalues tending 
to — oo. For a randomly chosen pair rj, u the kernel of the operator does not contain 
mean-free magnetic fields, and generically the only mean-free solution is b = 0. 

The processes bringing the system to a steady state thus can be viewed as adjustment 
of the flow to a conflguration allowing for a non-zero neutral mean-free magnetic mode. 
It is natural therefore to treat (1) as an equation in u. "Uncurling" it, one obtains 

TyVxb^uxb - TyVa, (3) 

where a is a scalar function (the constant factor rj is introduced for convenience). 

Consider separately the components of (3) parallel and perpendicular to b. Scalar 
multiplying (3) by b flnd 

(b- V)a = -b- (V X b). (4) 

The equation controls the magnitude of a magnetic fleld, whose direction is prescribed: 
Let is be a unit vector coUinear with b, then (4) implies 

(ij? ■ V)a . . 

^ "iB-(Vxi.)- 

If a magnetic force line is a closed loop (including the loops emerging due to spatial 
periodicity), then by virtue of (4) 



^ ifi • (V X b) = - ^{iB ■ V)a ds ^ 



(the parameter s on the curve is the distance along the curve from a flxed point on it), 
which can be also viewed as a constraint on the magnitude of the magnetic fleld (following 
from (5) ). 

The component of (3) perpendicular to b is accessed by cross-multiphcation of (3) by 
b, yielding 

r^b X (V X b + Va) = u|bp - b(u ■ b). 

The component of u parallel to b is not determined by (3), hence (u • b) remains an 
unidentifled arbitrary scalar function. We denote 

(u ■ b) ^ 
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implying 

u= (1 + 770;) b + Tye X (V X b + Va), (6) 

where 

e = b/|b|2. 

(6) and (4) together are equivalent to the equation for a neutral magnetic mode. The 
scalar field a satisfies the equation 

(b- V)q; + V- (e X (V X b + Va)) = 0, (7) 

equivalent to the solenoidality condition for the flow u. 

Now, in order to find the flow velocity (6), we need to determine Va, which we do 
employing the solenoidality condition for the flow. Generically and V x e are not 
parallel, and in this case (4) and (7) are equivalent to the equation 

Va = AiA + BiB + Cic, (8) 
where i^, ifi, ic is an orthonormal basis, 

. _ V X e - (b • (V X e)) is . _ . 

|V X e - (is ■ (V X e))iB| 

5 = -iB-(Vxb) = -|b|iB-(VxiB), (9) 

^_ (E/|b|)^-(b-V)a-V-(ex(Vxb)) 
|V X e- (is ■ (V X e))iB| 

We derive from (8) individual equations in A and C. 

The solvability condition for (8) is obtained by taking its curl: 

= VA X u + AV X u + VS X ifi + SV X ifi + VC x ic + CV x i^. (11) 

Scalar multiplying it by Ia, is and i^, one finds 

^ _ -CiA ■ (V X ic) - {iB ■ V)C - BiA ■ (V X ib) + (ic • V)^ 

iA • (V X i^) ' ^^^^ 

= Ais • (V X i^) + (ic ■V)A-^ + C'lB ■ (V X ic) - (i^ • V)C, (13) 

\o\ 

_ -Ale ■ (V X i^) + {iB ■ V)A - Bic ■ (V X \b) - {\a ■ V)B 
^ ~ ic • (V X ic) ' ^''^ 

where B is defined by (9). Substitution of (12) into (14) yields a second order differential 
equation along magnetic force lines, in principle, defining C. Initial conditions for this 
equations must be set on two-dimensional manifolds, transversal to magnetic force lines. 
They must assure geometric consistency: the solutions along closed force lines must be 
periodic. For force lines, intersecting with the boundary of the region occupied by the 
fluid, it is naturally to set the conditions on the boundary. The data can be provided on 
two manifolds, crossing a force line; in this case one obtains a boundary value problem 
for C. In turn, a can be found, in principle, from (10). This completes reconstruction 
of the flow. (The divergence of (8) yields an equation in a, which can be used to flnd 
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the potential itself.) Substituting A (12) and C into (13), one obtains an equation in b. 
Thus, not every solcnoidal field can be a magnetic neutral mode: The scalar consistency 
equation (13) constrains, together with the solenoidality condition, a neutral mode up to 
a scalar field. 

Implementation of this program can become particularly difficult in the presence of 
magnetic nulls, i.e., points, where magnetic field vanishes. (This is clear, of course, already 
from the definition of the vector field e, which becomes singular at the nulls). Topology 
of magnetic field with null points and its bifurcations during reconnections are studied in 
detail in solar magnetohydrodynamics [9, 10, 2, 17, 18] - they are presumed to be of funda- 
mental importance for occurrence of sudden explosive energy release events, solar flares, in 
the Sun's corona. In the vicinity of a null point magnetic field exhibits an approximately 
linear behavior controlled by the Jacobian ||96j/9a;j||. Solenoidality of the magnetic field 
implies, that the sum of the three eigenvalues of this matrix vanishes. Hence, generically 
it has two eigenvalues with real parts of the same sign, and an eigenvalue of an oppo- 
site sign. Consequently, one can identify a two-dimensional manifold of magnetic force 
lines behaving coherently - all approaching the null point or all departing from it (if the 
two eigenvalues have negative or positive real parts, respectively) and an one-dimensional 
manifold (a force line), exhibiting the behavior of the opposite kind. In the parlance of 
solar physics, the two-dimensional manifold is the fan, and the one-dimensional manifold 
the spine of the null (see Fig. 1 in [4]). Therefore, in our problem there are infinitely many 
characteristics (constituting the fan) , which must bring the same values of A and C to (or 
take the same values from) the null point, implying that the problem of consistency of the 
global solution for the fiow arises. The situation is further complicated by the fact that 
ifi is typically discontinuous at null points (its direction is not well-defined), and hence 
and are discontinuous as well. 

Thus, the presence of magnetic null points is likely to result in a discontinuity of 
the reconstructed fiow, but they are not the only source of troubles. More generally, our 
formalism becomes ill-defined at the points, where the magnetic field b is parallel to V x e. 
If a magnetic force line crosses the boundary at two points, a problem arises in satisfying 
the boundary conditions for the fiow at the two points. 

3. Axisymmetric magnetic neutral modes 

Equations (12)- (14) suggest that the complexity of the problem depends considerably on 
the geometry of magnetic force fines. For instance, reconstruction of the fiow is difficult, 
if force lines exhibit a chaotic spatial behavior. We consider here one of the simplest 
examples of an axisymmetric magnetic neutral mode 



(p, (p, z) being a cylindrical coordinate system and ip, i^, the respective unit vectors. 

Before we formulate the system of equations (12)-(14) in the variables A and C, which 
we need to solve in order to reconstruct the flow (6), we derive some useful properties of 
the basis i^, i^, ic- Curls of azimuthal and poloidal vector flelds independent of (p are, 
respectively, poloidal and azimuthal; this implies the orthogonality 



b = b{p, z)l 



Ia ■ (V X i^) = is ■ (V X is) = ■ (V x e) = ic ■ (V x ic) = 0. 



(15) 



By a simple calculation. 
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where ^ 



6' ^"'^ \VKy 

hence 

= X ic = hVn. 

Therefore, 

ic ■ (V X ia) = ic ■ (V/i X V/t) = 0, 

since none of the factors in the triple product has an azimuthal component. By vector 
algebra identities, 

U • (V X ic) - ic • (V X i^) = -V • (i^ x ic) = V • i^ = 0, 

implying 

U • (V X ic) = ic • (V X i^) = 0. (16) 

Now, scalar multiplying (11) by i^, i^ and ic and employing (15) (in particular, B — 0) 
and (16), one obtains equations 

= ^ifi • (V X i^) + (ic • V)A + C'lB ■ (V X ic) - {lA ■ V)C, (18) 

(which are now significantly simpler than (12)-(14) in the general case). Equations (17) 
and (19) are equivalent to 

C^C{p,z), (20) 
A = A{p,z). (21) 
For an axisymmetric magnetic field, (10) takes the form 

^ = -|V X elC - V • (e X (V X b)). 

Consequently, (20) and geometric consistency (27r-periodicity of a in (p) imply that 
a = a{p,z) is an arbitrary function (together with the relations (20), (21) and B — 0, 
this formally confirms a physically obvious fact, that a flow generating an axisymmetric 
magnetic field is necessarily axisymmetric), and (20) is superceded by 

V.(ex(Vxb)) 

|V X e| ^ ^ 

Now A must be determined from (18). We introduce characteristics {R{s),Z{s)) in 
the (p, z) half-plane; they satisfy the ODE's 

^^-h{R{s),Z{s))^{R{s),Z{s)), 
^^h{Ris),Zis))p^{R{s),Zis)). 
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Direct differentiation shows that the characteristics are isohnes of the scalar field b{p, z). 
Since along a characteristic 




(18) takes the form 



d_ 
ds 



{Ah) = /, 



where 



f = h{{iA-V)C-CiB-{Vxic)) 



and C is given by (22). Consequently, 



A{R{s),Z{s)) 



AjRiO), ZmhiRjO), ZjO)) + /q- fiRjs'), Z{s')) ds' 
h{R{s),Z{s)) 



(23) 



If a characteristic is a closed orbit of period S, geometric consistency implies that over 
this orbit 



Thus, we have determined Va and the flow (6) (to the extent this is permitted by the 
natural non-uniqucncss of solutions to (1) in u). 

The well-known Cowling antidynamo theorem states that generation of smooth ax- 
isymmetric magnetic fields (including steady ones) of finite total energy is impossible. Two 
proofs of the theorem (following [7] and [3]) are presented in [11]. The demonstrations 
rely on the equation of total magnetic energy balance derived for a smooth axisymmetric 
flow of incompressible fluid, provided the normal component of velocity vanishes on the 
boundary of the region where the fluid resides. To reconcile our results with the Cowling 
theorem, we note that the flow that we obtain will not satisfy some of these conditions. It 
may be singular on the circles, where b — 0, or k has extrema (and then e or h are singu- 
lar, respectively). If the volume occupied by the flow is bounded, it cannot be guaranteed 
that the normal component of the fluid velocity vanishes everywhere on the boundary 
(or, alternatively, enforcing this condition creates a discontinuity in the flow). Hence, 
the standard procedure employed to establish the total magnetic energy balance equation 
will reveal additional sources of magnetic energy, which emerge because the flow is not 
smooth or the surface integral representing the contribution of the advective term does 
not vanish; under such circumstances the Cowling theorem is unapplicable. 

We have presented the analysis of this section mainly as an illustration of how the 
proposed formalism might be applied to reconstruct flows for less trivial magnetic field 
configurations. However, in addition, it provides useful information in regard to the 
following technical issue: Although we have stated at the end of the previous section that 

(13) is a constraint for a neutral magnetic mode, we have not yet produced any evidence, 
that the three equations (12)-(14) arc independent. Eqns. (17)-(19), which we have derived 
considering this particular example, demonstrate that (13) is not a consequence of (12), 

(14) and solenoidahty of magnetic field. 




(24) 
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4. Concluding remcirks 

We have shown in Section 2 that reconstruction of an incompressible flow (6) from the 

structure of a magnetic mode consists of solution of equations (12) and (14) in A and 
C, followed by solution of (10) in a. These equations are ordinary differential equations 
along magnetic force lines; thus, the problem becomes complex, if the force lines exhibit 
a chaotic behavior. For a solenoidal vector field to be a neutral magnetic mode, it must 
satisfy the constraint (13). 

Substitution of (6) into the momentum equation 

i/V^u + ux(Vxu)-bx(Vxb)-Vp + F = 

yields an equation in b: 

i/V^((l + 770;) b + 776 X (V X b + Va)) 

+ r){{l + r]a)h + r]e x (V x b + Va)) x (V x (ab + e x (V x b + Va))) 

+ r]{ah + e X (V X b + Va)) x (V x b) - Vp + F = (25) 

comprising a closed system of equations together with the solenoidality condition (2). 
Relation (13) now becomes a constraint on the acceptable fluid forcing F. 

Analysis of the dependence of steady or evolving magnetohydrodynamic systems on 
small viscosity and magnetic diffusivity is a notoriously difficult problem. The structure 
of (25) may turn out to be advantageous for the study of asymptotics of MHD steady 
states, when the force F is of the order of small quantities ~ 77, as it is in nonlinear 
dynamos with energy equipartition [5, 6]. (The form of the scalar factor in front of b 
in (6) has been chosen so that all terms in (25) were in this case of the same order of 
smallness.) 

In Section 3 we have considered an example of the reconstruction problem for axisym- 
metric neutral magnetic modes. This particular case has proved to be highly degenerate: 
the denominators in (12) and (14) vanish identically, and the respective components of 
(11) just testify that Va is an axisymmetric vector field. Relation (13) does not constrain 
further the structure of the magnetic field, but rather defines, by (23), the component A 
of Va. Initial conditions A{R{0), Z{0)) for solutions (23) of (18) along characteristics can 
be chosen on curves in the (p, z) half-plane, which are transversal to magnetic force lines. 
The azimuthal component of the flow velocity, (1 + r]a)b, is an arbitrary axisymmetric 
scalar field (in this case it is controlled neither by the magnetic induction equation, nor, 
due to independence of (f, by the solenoidality condition). Thus the reconstructed flow is 
unique up to the data which must be specified on two-dimensional manifold(s) (the scalar 
field a on the (p, z) half-plane) and on one-dimensional curve(s) on this half-plane (the 
initial conditions A{R{0), Z{0)) ). 

The initial data must be smooth so that the resultant field A had no singularities. If the 
topology of isolines of the magnitude of magnetic field b is non-trivial, the smoothness of 
the initial data is insufficient; for instance, geometric consistency requires that the integral 
(24) over any closed magnetic force line vanishes. If the axis of symmetry intersects with 
the volume occupied by the fiuid, axisymmetry gives rise to another problem: regularity 
of the magnetic field implies 6(0, z) = 0; consequently, the term e x (V x b) in (6) tends 
to infinity for p ^ 0. Thus, the flow is non-singular only, if initial conditions for A 
compensate for this singularity. 
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